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ABSTRACT 


A Green's  function  method  has  been  used  to  obtain  an  expression 


for  the  mean  free  path  of  a Rayleigh  wave  propagating  along  a planar 
free  surface  of  an  isotropic  elastic  continuum  and  scattered  by  a 
mass  defect.  The  change  in  density  associated  with  the  mass 


defect  is  assumed  to  be  Am6(x-x  ) where  x is  the  position  vector 


of  the  defect  and  Am  is  the  mass  change.  The  Green's  function  is 
evaluated  for  an  isotropic  elastic  continuum  with  a stress-free 
planar  surface.  Using  this  Green's  function,  the  continuum  equations 
of  motion  are  formally  solved  for  the  particle  displacement  of  the 


Incident  wave.  The  Poynting  vectors  are  then  calculated  for  the 


incident  wave  and  the  scattered  wave 


scattered  wave  Poynting  vector  are  obtained  in  the  asymptotic 
limit  of  large  distance  from  the  mass  defect.  The  mean  free  path 
is  then  obtained  from  the  ratio  of  the  magnitudes  of  the  incident 
Poynting  vector  and  the  asymptotic  scattered  Poynting  vector.  The 
results  are  compared  with  those  of  other  workers. 
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I.  Introduction 


The  development  of  electronic  devices  using  surface  elastic 
waves  has  stimulated  considerable  Interest  In  the  fundamental 
properties  of  these  waves.  Important  among  these  properties  Is 
the  rate  of  scattering  of  surface  waves  due  to  their  Interaction 
with  various  types  of  defects  and  with  other  phonons.  The  scatter- 
ing rate  can  have  a significant  Influence  on  the  performance  of 
devices. 

Theoretical  Investigations  of  the  anharmonlc  scattering  of 

Rayleigh  waves  have  been  carried  out  by  a number  of  workers. 

Using  a Green's  function  procedure,  Maradudln  and  Mllls^1^  treated 

the  Isotropic  case  using  a lattice  dynamical  model  and  found  that 

4 

the  damping  constant  Is  proportional  to  iurT  at  low  temperatures, 
where  u)R  Is  the  Rayleigh  wave  frequency  and  T Is  the  absolute 
temperature.  These  results  were  extended  to  anisotropic  crystals 


by  King  and  Sheard. 


On  the  experimental  side  the  idrT  depen- 


dence was  verified  by  Salzmann,  Pllenlnger,  and  Dransfeldv  ’ for 
several  surfaces  of  quartz. 

The  effect  of  defects  on  the  damping  of  Rayleigh  waves  was 

(4 ) 

Investigated  theoretically  by  Steg  and  Klemens  for  the  case  of 
a point-mass  defect  having  a mass  change  Am.  They  used  perturba- 
tion theory  and  found  that  the  relaxation  rate  Is  proportional  to 

(Am)  u)R  . Somewhat  later,  Sakuma'  ’ ’ re-examined  the  problem  

using  the  Chew-Low  scattering  formalism  and  the  complete  set  of  < 

SnMm 

normal  modes  for  a semi- Inf lnlte  Isotropic  elastic  medium  con- 

(7)  5 ~ ~~ 

structed  by  Ezawa.  Sakuma  confirmed  the  u)R  dependence  of  the 


lOTiBuiir.'M!'.!?!';!'  iris  ; 

AVAil  a.iO  i . iAL  1 


□ □ 
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scattering  rate  at  sufficiently  low  frequencies,  but  in  addition, 

found  resonance  structure  for  a defect  with  a lighter  mass  than 

the  host  atoms.  This  work  has  been  extended  to  the  case  of  ran- 

ts) 

dom  density  fluctuations  on  the  surface  by  Nakayama  and  Sakuma. 

In  the  present  paper  the  Green's  functions  for  a semi- 
infinite  isotropic  elastic  continuum  with  a stress-free  planar 

(9) 

surface ' 7 are  used  to  calculate  the  inverse  attenuation  length 
of  a Rayleigh  wave  scattered  by  point-mass  defects.  In  Section  II 
formal  expressions  are  derived  for  the  amplitude  of  the  scattered 
wave.  In  Section  III  these  expressions  are  evaluated  explicitly 
in  the  asymptotic  limit  of  large  distance  from  the  mass  defect. 

In  Section  IV  the  inverse  attenuation  length  is  evaluated.  Nu- 
merical results  are  presented  in  Section  V.  A discussion  and 
comparison  with  previous  work  are  given  in  Section  VI.  As  we 
shall  see,  our  results  concerning  the  resonance  structure  differ 
significantly  from  those  of  Sakuma. 
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II.  Scattering  of  a Rayleigh  Wave  by  a Mass  Defect 


In  the  presence  of  a mass  defect  situated  at  the  point 
xQ  - (0,0,x0g)  In  a semi- Inf lnlte , Isotropic,  elastic  medium  occu- 
pying the  half-space  Xg>0,  the  equations  of  motion  of  the  medium 
can  be  written  In  the  form 

_p  .ii  u + e iik  + e c * y 

^ dxp  axv  puv  a^v 


P S L^(x,t)  ^ - im  6(x  - xQ) 


(2.1) 


where  utf(x,t)  Is  the  or  Cartesian  component  of  the  displacement 
field  at  the  point  x at  the  time  t , p Is  the  mass  density  of  the 
medium,  and  Am  is  the  Increase  in  the  mass  of  the  medium  due  to 
the  Introduction  of  the  defect.  In  writing  Eq.  (1)  we  have 
assumed  that  the  elastic  moduli  v(x)j  are  position  dependent 


and  are  given  by 


C«0uv  (x)  ” 0(x3)  C«0uv 


(2.2) 


where  the  {caguV}  are  the  usual,  position- Independent  elastic 
moduli  of  the  medium,  and  0(Xg)  is  the  Heaviside  unit  step  function. 
We  now  Introduce  a Green's  function  G _(x,x;t-t')  as  the 

®p 

solution  of  the  equation 

Ti  L (x,t)  G . (x  ,x# ; t - t')  - 6 6 (x  - x')  6 (t  - t')  , (2.3) 

U ®U  UP  ®P 

subject  to  outgoing  wave  or  exponentially  decaying  wave  conditions 
as  Xj  4 a, 
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In  terms  of  this  function  we  can  rewrite  Eq.  (2.1)  as  an 
Integral  equation: 

ua(x,t)  - u^o)(x,t)  + ^ S Jd3x'Jdt'G  (x,x#  ;t  - t')  x 


X 6(x'  - X ) -n  U (x'  ,t')  , 

° at'  p 


(2.4) 


where  u^°^(x,t)  Is  a solution  of  the  corresponding  homogeneous 
equation, 


5 Lau^,t)u^o)(x,t)  - 0 


(2.5) 


U “U 

and  In  the  present  context  represents  a Rayleigh  wave  propa- 
gating along  the  surface  Xj  - 0 of  the  semi- Inf lnlte  elastic 


medium. 


With  the  Fourier  decompositions 


u_  (x , t ) - u (x  ,«j)e” liwt  , u<o)(x,t)  - u^o)  (x  ,ui)e” ltt,t  (2. 


6) 


G«p^x,x  ;t~ t “ J "Stt  Gag(x’x';n>e 

,2 


(2.7) 


$ Gap(5^';n)e-10(t“t,) 

OaB(i,J';n)  - J eU«  ■ (x>  ■ *<)goB(ii1n|x3Xp  (2.8) 


where  X|  - (Xj.Xg.O)  and  £||  - (k^kg.O),  we  can  rewrite  Eq.  (4) 


as 


ua(x,«)  - u£o)(x,u>)  - ^ Gflrg(x,xo;uo)ug(xo,UJ)  . 


(2.9) 
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To  solve  Eq.  (2.9)  we  set  x - xQ  and  use  the  fact,  established 
in  the  Appendix,  that  ls  » diagonal  matrix.  In  this 

way  we  obtain  the  equation 

uor<V“l)  * “«°)<V“|)  ■ <210) 


with  the  solution 


ua(x0,u>)  - 


1 + Gaa(*o’*o;Uj) 

o or  or  o o 


(2.11) 


When  we  substitute  this  result  into  the  right  hand  side  of  Eq. 
(2.9),  we  find  that  the  amplitude  of  the  scattered  displacement 
field  ls  given  by 

uis)(x,W)  s u (x  ,u> ) - u^o)(x,tu) 

Of  « or 


bm 


ft  Go,s(J'*o;“)  VV“)  “v°)<;o’*>  • <212) 


where  the  scattering  matrix  is  given  by 


Vv*)  * 


1 + 


— ' . 6;p — 

w!G  (i  ; W"0 

p 'xo  ,xo  ' 


(2.13) 


For  the  system  under  consideration  it  is  also  the  case  that  the 

nonzero  elements  of  this  matrix  obey  the  relations  T^(x0,i»))  ” 

T„„(x  ,u>)  * T00  (x,.U) ) . With  the  aid  of  Eq.  (2.8),  the  scattered 
a o jj  o 

field  takes  the  form 


- - am!  r d2k„elS«  • 5<  s.^CV1-3-03)  uB(°>(; ..)  . 
* 4rr%  J P Dp(xo;w)  * ° 


(2.14) 


The  Fourier  coefficients  -[g^  (It||  uu  | x3x3 )}  can  be  expressed  in 
terns  of  another  set  of  coefficients  {dap(klialx3x3)}  by 

t«0  (^»  * 1*3X3)  - S <^v(ktttu  1x3X3)  Sua(fc,i ) Svp(fc„)  , (2.15) 


where  the  real,  orthogonal  matrix  S(&||)  is  given  by 


*i  K 0 


sd„)  - -fc2 


(2.16) 


The  nonzero  elements  of  the  tensor  dap  (k||  t»  | X3X3  ) , viz.  d11,d13 

d02 , d31,  d33  have  been  calculated  recently  for  a semi** Infinite, 

(9) 

isotropic,  elastic  medium  occupying  the  upper  half-space  x3 > 0. 
Consequently,  all  of  the  functions  entering  the  right  hand  side  of 
Eq.  (2.14)  are  known,  and  the  scattered  field  can  thus  be  calcu- 
lated. We  will  return  to  this  aspect  of  the  problem  below. 

For  the  Incident  wave  we  assume  a Rayleigh  wave  propagating 
in  the  positive  Xj-direction.  If  we  write  the  amplitude  u^°^(x,m) 
in  the  form  . , . 


ui0)(x,u>)  - uf°^  (It|f0*u)|x„)e 


i£«(o)  • Xu 


(2.17) 


with  k|f°^  - (k^°\  0,  0),  the  amplitudes  (S|f°^u> |x«) 


are  found 


to  be 


»<°>  <*,<0W3>  - A{8'B^3 -(l-lS-j)  e'P‘‘3}  (2.18&) 


“p*x. 


* 
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Uj0*  (<|(o)«.|x3) 
u3o)  (k|fo)u)|x3) 

where  A is  an  arbitrary  amplitude,  and  are  the  speeds  of 
longitudinal  and  transverse  waves  in  the  isotropic  elastic  medium, 
respectively*  and  cR  is  the  speed  of  Rayleigh  waves,  which  is 
obtained  from  the  equation 


- 0 


- t(l-A)  *{® 


(2.18b) 

-*tx3  , 

-2 g—  ]•  , (2.18c) 

x 2c 


The  other  quantities  entering  Eq.  (2.18)  are  defined  by 

2 , . 2 i 


U) 


- 6,  - »<•>(!  - . (2.20) 


The  results  expressed  by  Eqs.  (2.17)  and  (2.18)  have  the 
consequence  that 

(s),-*  . hw2  ui0)  <k|f°)t*,lxo3)  .(1),-*  , „ 

u ' (x ,u> ) “ - (x;x03|uj)  - 

4tt  p D1(xo;t«) 

. 2 ui°^  (it|f°^U)  | xn„) 

_ tm l_  _? ! L23_  U) 


4tt  p D3(xq;uj) 


(2.21) 


where 


1^*  (x  ; x03|uj)  - Jd2k||elk|1  * X"  $11  u>lx3xo3>  • (2.22) 


The  resonance  denominators  Da(x0  ; uj)  are  evaluated  in  the  Appen- 
dix, so  that  it  is  only  with  the  integrals  1^  (x  ; x03|uj)  that  we 
will  be  concerned  in  what  follows. 


7«r*K. 


— ’• 
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III.  The  Scattered  Displacement  Field 


In  this  section  we  obtain  the  asymptotic  behavior  of  the 
integrals  1^®^ (3  ; xo3|uj)  defined  by  Eq.  (2.22),  for  3 far  from  the 
point  (0,0,  xQ3) . In  the  following  section  these  results  will 
be  used  to  obtain  the  elastic  Poynting  vector  of  the  scattered 
displacement  field,  and  from  the  latter  the  attenuation  length 
of  Rayleigh  waves  due  to  scattering  by  mass  defects. 

Since  the  results  of  this  section  are  central  to  the  calcu- 
lations ih  the  rest  of  this  paper,  we  present  their  derivation  in 
some  detail. 

We  begin  by  substituting  Eq.  (2.15)  into  Eq.  (2.22): 
l£S)<*  i *03 1**  - 2 K*.*1*1  ' 

(3.1) 

When  the  explicit  expressions  for  the  matrix  elements 

are  employed  in  Eq.  (3.1),  the  six  integrals  1^(3  ; xQ3|uu)  which 

are  required  for  the  determination  of  the  scattered  wave  are 


given  by 

l[l}a  ; x03 1***)  - Ii  + *2  (3.2a) 

I^X)(3  ; x03|<h)  - l4  - l5  (3.2b) 

; x03|(tt)  - I?  (3.2c) 

j[3)(3  ; Xq3  |im ) - I3  (3. 2d) 

l£3)(x  ; xQ3|m)  - Ig  ' (3.2e) 
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x;Xo3|uj)  - Ig, 

(3.2f) 

simpler,  auxiliary  integrals  1^,  ..., 

Ig  are  defined  by 

J d^iie1^"  '*l,cos2cpdj1(k||U)|x3x03) 

(3.3a) 

J d^ye1^11 ’X||sin2cpd22(kM(u|x3x03) 

(3.3b) 

J d^ye1^11  *x,lcoscpd13(k|uu|x3x03) 

(3.3c) 

J d2kyeik|t  ’^“coscpsincpdj^kinjuIxgXQg) 

(3.3d) 

j*  d2k||eik|1  x,,coscpsincpd22  (kyuj  |x3xQ3) 

(3.3e) 

J d^kye1^11  *x,lsinqpd13(k|ou  |x3xQ3) 

(3 . 3f ) 

J d^ky  e1^11  *X|lcos(pd31(kMu)|x3x03) 

(3.3g) 

J d2l[||elk''*“  d33(k,ul|x3x03)  . 

(3.3h) 

Ve  next  represent  the  vector  x t,  as 
X||  - x ||  (coscpg , sincps,  0) 


(3.4) 


and  evaluate  the  Integrals  over  cp  In  Eq.  (3.3)  in  the  limit  of 
large  x(|  by  the  method  of  stationary  phase.  The  results  are 


h ~ d?)*  e_i  4 coa\  Jn 
l2  ~ d?)  e * sin2«Ps  j22 


x3  ~ dr) 


i -ij 


e COSm  J 


^ J13 


(3.5a) 

(3.5b) 

(3.5c) 


~T"  HM 


i -l5 


~ (If)  9 s incPgCos cps  Ju 

i -i  n 

X5  - dj)  e ? sincpscoscpg  J22 

i -i  n 

h ~ d?)  e * sln®s  J13 


J -i  n 

J7  ~ df)  e * coacp3  J31 


*s  ~ df) 


i -i*4 

e 


'33 


where 


CD 

Ja0(*;XO3  " J dkk^e  lldaB^ka,lx3X03^* 
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(3 . 5d) 

(3 . 5e) 

(3.5f ) 
(3. 5g) 
(3.5h) 

(3.6) 


In  what  follows  we  focus  our  attention  on  the  integrals  ^ (x ; 3Cq3  | cu > 

The  simplest  of  these  is  J 22 , and  its  evaluation  illustrates 

many  of  the  ideas  to  be  used  in  the  evaluation  of  the  remaining 

integrals.  From  the  results  of  Ref.  9 we  have  that 


i 


J22(x;x03lu,)  " ‘ -72  J tr 

ZC  . t 


ikxii 


x e 


D 


t o 

— Of  * ( x0  • X. 


t *3  03 


| "at^X3+X03^ 


+ e 


where 


(3.7) 


(3.8) 


The  correct  analytic  continuation  of  from  the  region  with  k > U)/ct 
to  the  region  k<ou/ct  is  achieved  by  taking  the  branch  cut  al'ong 
the  negative  real  axis,  and  assuming  to  to  have  an  infinitesimal, 
positive  imaginary  part. 
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Each  of  the  two  integrals  in  Eq.  (3.7)  is  of  the  form 

A - ? dk  — ef(k) 

_ “t 


where 


f (k)  - ikxii 


■ $ 


(3.9) 


(3.10) 


with  R either  |x3  - x3 | or  (x3  + xQ3).  We  expand  f(k)  about  its 


stationary  point  kQ  , 


f(k)  - f(kn)  + \ (k-k0)V(k  ) + ••• 


where 


k - JL 


o c 


t (xj  + R2) 


f(ko)  " 1 ^ (x?  + R2)‘ 
2 2 

C t (xJ+R2) 

£'<V  -15-71 


(3.11a) 


(3.11b) 


(3.11c) 


(3. lid) 


Since  kQ  is  in  the  interval  (o,  «)  (in  fact*,  it  is  in  the  interval 
(o,  u>/ct)),  we  can  replace  the  integration  interval  (o,«  ) by  the 
infinite  interval  ( - »,  ®) , in  the  limit  of  large  X|(  and  x3,  with 
an  error  which  is  of  higher  order  than  the  terms  we  retain.  In  this 


way  we  obtain 


i f-  (xjf  + R2)  . 

e * J dk  e 


n O 3/2  n 
iCt(xjf  + R2)  k2 

"5  uj  12 


(2nx„)i  -i  | 1 7%  (x»  +R  > 


(*?  + R2) 


(3.12) 
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We  now  make  the  assumption  that  xQ3 
2 2 ^ 


with  x - (X||  + x3) 
expanded  in  powers  of  (xQg/x)  with  the  result  that 


is  small  In  comparison 

£ 


2 2 - 

The  two  values  of  (xj  + R ) can  therefore  be 


2 2 ^ x3  t xll  o 

(x,|  + R ) - x ± — x03  + i -J  x03 


+ • 


(3.13) 


where  the  upper  (lower)  sign  obtains  when  R - (Xg  + xQ3 ) (R  - |xg  - xQ3 
We  will  retain  the  term  linear  in  x^g  in  the  exponential  factor  in 
Eq.  (3.12),  but  only  the  term  of  zeroth  order  in  xQg  in  the  deno- 
minator of  the  prefactor.  The  justification  for  this  assumption, 
which  it  should  be  emphasized  is  convenient  but  not  essential,  is 
the  following.  In  the  expressions  (2.21)  for  the  amplitudes  of  the 
scattered  displacement  field  the  integrals  (x;xQg  |uj)  appear 

multiplied  by  the  amplitude  u^°^  (ic|f°^uu  |xQg)  of  the  incident  Rayleigh 
wave  at  the  impurity  site.  From  Eqs.  (2.18)  we  see  that  these  latter 
amplitudes  decay  exponentially  with  increasing  xQ3,  so  that  it  is 
only  for  xQ3  p”1  , where  Pt(<  P^)  is  defined  by  Eq.  (2.20),  that 


* 

any  significant  scattering  of  the  incident  Rayleigh  wave  by  a mass 
defect  can  occur.  In  evaluating  the  scattered  field  far  from  the 
impurity  site,  we  assume  that  kQx  » 1,  where  kQ  is  the  wave  vector 
of  the  incident  Rayleigh  wave  and  is  comparable  to  the  wave  vectors 
of  the  scattered  waves,  as  we  will  see.  If  kQx  » 1 and  xQ3  p”1, 
it  follows  that  xQ3  « x,  and  the  .approximations  we  are  making  should 
lead  to  little  error  in  the  scattered  displacement  field.  To  this 
approximation  we  obtain  finally  for  the  integral  J22(x;xQg  |uu)  the 
result  that 


J22(x;xQ3  | a; ) ~ - 


03 


± (2  tt  x, ) 

"2 x 

ct 


* -i?  ‘f1 

- e * e * 


UU 


008  T2-T  *03-  <3-14> 


). 


iirvist 
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It  follows,  therefore,  that  the  Integrals  I2  and  I5  are  given 


asymptotically  by 


t 2tt  e 

l2 2 sin  CPS  ” 

ct 


iktx 

cps  S-Y-  cos(ktx03cosflg) 


ik.  x 
t 


where  we  have  used  the  fact  that 


X||  - x sin0  , x«  • x cos  0S 


and  have  defined 


k - 

kt  ct  ‘ 


(3.15a) 


Ij ^ siacpgCoscpg  — * — cos ^ktx03cos®s^  * (3.15b) 


(3.16) 


(3.17) 


We  now  turn  to  the  remaining  integrals  ^ (x ; Xq3  | uj  ) with  af,p-l,3. 
The  Green's  function  dag(kcu |x3xQ3)  for  - 1,3  can  be  written  as^9^ 


/T  \ “®  • | | 

dae(ku'lx3x03)  ” dap  ^ku,^e 


_a-tX3 


-ot+x. 


+ d^o  (kou)e 


“®t  X3“X03 


"“tX3 

e t 3 f„(5) 


Kgto)e 


““-t,x03 


+ d<f  > (k«)e"°'tX°3  } 


where  c*t  has  been  defined  in  Eq.  (3.8),  while 

/ 


•*  • (-'  ■ 5) 


(3.18) 


(3.19) 


r±  * 


4«tcr  c^k2  ± (u,2  - 2c2k2)(a2  + k2) 


4atcr^(ou  -2c^k  ) 


(3.20) 
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We  have  factored  out  the  coefficient  r“  explicitly  in  Eq.  (3.18), 
because  r”1  has  a simple  pole  at  k ■ kR  ■ (uu  + io)/cR,  where  cR(<  ct) 
is  the  speed  of  Rayleigh  surface  waves.  Explicit  expressions  for 
the  coefficient  functions  d^(kuj)  i - 1,2,..., 6 are  given  in  Table  I, 

The  asymptotic  forms  of  the  two  integrals  containing  d^P(kou) 

\ «P 

( 2 \ 

and  dv  '(kcu)  can  be  obtained  immediately  with  the  aid  of  the  results 
we  have  already  obtained.  In  the  former  case  it  is  necessary  only 
to  replace  ct  by  c^  in  the  preceding  analysis.  Thus  we  obtain 
immediately  that  the  integral  ^ 

J^xjXQglw)  s J dkk*  e^'d^  (k<u)e  (3.21a) 


with  (or,p  - 1,3)  is  given  asymptotically  by 


.m(j  , , 1 -i  f e^VlV03  ®. 

J<*(3  <x'x03l“,)  ~ T -T2 e x 

2cl 


sin2  0. 


where 


sin  0 cos  6 
s s 


sin0„cos0_  cos  0, 


k ■ JL 
* 


(3.21b) 


(3.22) 


In  a similar  fashion  we  find  that  the  integral 


\ 

J<f(5:Xo3|»)  > jdk  d«>(k»)e'"t|X3'X°3 


(3.23a) 


is  given  asymptotically  by 


, (2tTX||  )2  -i  -r  . 

JV<x'*03l">  * 5 


-iktx03coses 


V 


2 

cos  0, 


-sin  0 cos  0 
s s 


-sin0  cos0_  sin  0„ 

s s s 


(3.23b) 


We  turn  now  to  the  integrals  which  contain  the  factor  r~ 
in  their  Integrands.  We  will  work  out  one  of  them  explicitly,  and 
will  simply  quote  the  results  for  the  remaining  three.  The  integral 
we  consider  is 

**03  I-)  - j * «lk“  ' V3  £ o'3.'  V°3  • (3.24, 


To  evaluate  this  integral  in  the  limit  of  large  x we  divide  the 
range  of  integration  (o,  ®)  into  two  parts:  (o,  uu/c^)  and  (uu/c^,  ®) , 
and  consider  the  resulting  two  integrals  separately. 

A.  The  Contribution  from  the  Interval  (o,  uu/c^). 


We  evaluate  the  integral 

(ju/c 


t (3a)  .j  I n _ r Hlr  i 3 k2  .(3) 

Ja0  J ^ e r dap 


f d~a  (kw)e 
r+  ap 


~alX03 


(3.25) 


by  the  method  of  stationary  phase.  If  we  denote  by  kQ  the  stationary 
point  of  the  function 


f (k)  - ikxi 


we  find  that 


-(*2-4)  X3  “ f(ko>  + * f'<V(k- 


ko>  + 


k - J*L_1  < JL 
° c*  x c* 


f(kc)  - iktx 

c,  x3 

f'<ko>  -*1¥7 

3 


(3.26) 


(3.27a) 


(3.27b) 


(3.27c) 


We  therefore  obtain  for  the  large  x limit  of  the  integral  (3.25) 
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J^a)(2;x03|u)) 


-[fi  J 

“ GO 


1 °l  *3  k2 

2 ~ “5  k 


>*  -i? 


ik.x 


(2tnc| )2  -i  J e“t*  ik^xO3cos08 


(1-  2X2sin20,a)2-  4X3sin20  cose  (1-  x2sin29_)* 

x 2 —t  x 

(1- 2\2sin20e)  + 4x3sin20  cos0  (1  - X2sin20„) 

S S3  O 


-sin2  0. 


-sin0  cos0„ 
s s 


sin@„cos  0„ 
s s 


cos20_ 


(3.28) 


where 


X - ct/ct  < 1 . 


(3.29) 


We  have  used  the  results- of  Table  I in  obtaining  this  result. 

B.  The  Contribution  from  the  Interval  (ou/c^,  ®) 

To  evaluate  the  contribution  to  (x;xQ3  |uu)  from  the  integra- 

tion range  (m/c^)  < k < ®,  we  regard  k as  a complex  variable  and 
consider  the  contour  integral 


?aB  ^x»x03lu^  “ J dke 


ikxj|  — a.x«  . J /o\  — ®/x03 

4 <3-30) 

r+  orp 


where  the  contour  C.  is  shown  in  Fig.,  l.  The  integrand  can  be 
made  single  valued  inside  and  on  this  contour  by  a proper  choice 
of  branch  cuts.  Then  by  employing  the  same  arguments  that  were 


used  in  Ref.  9 in  connection  with  the  evaluation  of  a similar 


r 
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Integral,  In  the  limit  as  the  circular  portion  of  the  contour 


recedes  to  infinity  the  integral 


J^b)(x;x03|u))  - J dk  e 

i u/c , 


ikX||  - «.X3  ki 


**3  S!  d<3>(k«,)e’V°3 
r+  erp 


(3.31) 


is  given  by  2 tt  i times  the  residue  at  the  simple  pole  the  integrand 

possesses  at  the  zero  of  r+  at  k - kR  ■ (uu  + io)/cR.  In  the  vicinity 

(9) 

of  this  pole  we  have  the  expansion v ' 


f-  - - §5  ;k  - kry  + 


(3.32a) 


where 


(-&)M  • J "J  1 *-;j 


32b) 


Consequently,  we  obtain  immediately  that 


'Sb,<s**osl->  (£)  e 


5 ikRxll  -et<*3+x03) 


Hl’-HlVl)1  ' -H) 

(3T33) 

In  the  same  way  we  obtain  the  asymptotic  behaviors  of  the 
remaining  integrals.  We  omit  the  details  and  merely  present  the 
results.  We  have  that 
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- J * ."v>  # -“W* 


-Qf^X 


03 


l (2ttxw) 
2c 


i 


~ -j g cos0asin0ge 


-t  J .V  1 ^ x03(1-  ^3il>2V 


4Xcos0a(l  - 2\2sin20a) (1  - X2sln20a) 


i 


(1  - 2X2sin20s)2  + 4X3sin20acos0a(l  - x2ain20#) 


tan0. 


\ 


Xsin20. 


cos0_(l-  X2sin20„) 

9 9 


T 


Xsin0, 


(1-  X2sin20a) 


T 


21 


4€(e  )cose  |\2-  ain2e  |'(1-  2sin*e_) 


a s 


4€  (0  )sin20  cos0  I X2  - sin2 9 1 + (1-  2sin20  )cos2  0„ 

9 S 9 1 ' S S 


sin  9, 


€ ( 0 ) | X2  - sin20  | 


sin  0. 


- tan0. 


€(6a)|x2-»in2e.| 


tan  0, 


2tt  i / uu  \5/2  ikRx"  " B-tx( 


2tt  i ( uj  \w/  R 11 

+ — s e e 

R<i)  VCR7 


HJHl 


M)  -H)H) 

H)‘ 


where  we  have  introduced 


€(0a)  - 1 0 <sin0g  < X 


- i X <ain0g  <1  ; 


o + 

j (2tt*h  -if9lktx  1 F;  *O3ooa0s 


‘^*03 


(3.35) 


(3.36) 


a 


4€(0g)sin20gcos0ajx2- sin20g|  - (X  - 2sin2 0g)cos20£ 
4€(0a)sin20scos0a|\2- sin20g|  + (X  - 2sin20g)cos2  0s 


o 

- cos  9. 


sin 8 cos 0„ 
s s 


- sin0_cos0_ 

S S 


sin2  0, 


(#-) 
2Ruj  'CR7 


5/2  ikRx„  -et<x3  + x03) 


.(3.37) 


"H) 

MHJH)  4-S*  ■ 

\ t 


The  results  of  this  section  when  substituted  into  Eq.  (2.2X) 
enable  us  to  write  the  amplitude  of  the  scattered  displacement 
field  u(s)(x,u,)  as  the  sum  of  four  contributions: 

u(s)(x,uu)  - u^(x,u>)  + u<tp>(x,u,)  + u^ts^(x,uu)  + u(R)(x>UJ)  . 

(3.38) 

In  Eq.  (3.38)  u^(x,uu)  is  the  amplitude  describing  the  scattering 
of  the  incident  Rayleigh  wave  into  bulk  longitudinal  waves; 
{T(tp)(x,u,)(iT(ts)(x,u)))  is  the  amplitude  for  scattering  into  bulk 
transverse  waves  of  p-polarization  (s-polarization)  ; and  vt^R\i,uj) 
Is  the  amplitude  for  scattering  into  other  Rayleigh  waves.  The  ex- 
plicit expressions  for  each  of  these  amplitudes  have  the  following 
simple  forms. 
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-*(/.)  x x Anui  e 

a)  uvc;(x,m)  - - 2 — 

4 TT  pc^ 


ik^X 


(3.39a) 


u.  - 


ul°^  lx03 ^ 


4 D1(*0;u“) 


sin  0gcoscpgF  ^ ( 0g , XQ3 ) + 


u30>  (^|[0>u)  |x03) 

d3(x0;«>) 


cos0sF2^ 0s,xO3^  » 


(3.39b) 


where 


Fl(es*x03)  “ e 


“iktxO3cos0s 


iktxO3cos0s  (1“2\2sin2eg)  -4X3sin2egcos0g(l-X2sin29g) 


2 o 

( 1-2  X2s  in2©^ ) +4X3s  in2  0gcos  0g  ( 1-  X2s  in2  0S ) 


iktxO3(1-x2sin20s) 


4 Xcos 0 ( 1- 2 X2s in2  0g ) ( 1- X2s in2  0g ) 


( 1-2  X2s in2  0g ) +4x3s in2  0gcos  0g ( 1- X2s in2  0g ) 


(3.39c) 


. /fl  , -iktxO3coa0s 

F2(0s*xO3)  " e + 


2 1 
ik£x03c°s eg  (l-2x2sin20g)  -4x3sin29scos0g(l-X2sin26s) 

» 5 ; 

( 1-  2 X2s  in2  0O ) +4  X3  s in2  9gcos  9g  ( 1-  X2s  in  2 0Q ) ' 


iktx03(1"x2sin2es) 


4X2sin2 0g(l-2X2sin2 0g) 


(l-2x2sin20g)  +4x3sin29gcos  0g(l-X2sin20g) 


(3 . 39d) 


ik^x 


b)  U(tp)(x,aj)  - -A  u<p)  -2— — [x^os 9scoscps  + 

4tt  p 

+ XgCos  9gsinqpg  - x3sin0g  J , 


(3.40a) 


I 
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(3.42b) 


We  now  turn  to  the  determination  of  the  elastic  Poynting 
vector  of  the  scattered  displacement  field,  and  from  it  tho  in~ 
verse  attenuation  length  of  a Rayleigh  wave  in  the  presence  of 
point  defects. 
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IV.  Evaluation  of  the  Inverse  Attenuation  Length 

The  first  step  in  the  evaluation  of  the  inverse  attenuation 
length  is  the  calculation  of  the  elastic  Poynting  vector  for  the 
incident  and  scattered  waves.  In  their  paper  concerning  the 
attenuation  of  Rayleigh  waves  by  surface  roughness,  Maradudin 
and  Millsv  derive  the  following  expression  for  the  complex 
Poynting  vector  Qc  of  an  elastic  wave  propagating  in  an  isotropic 


medium, 


c - (c? > d,  cS> 


(4.1) 


where 


l - - ip*l{oJ(7-3)  - - 


<4-2) 


c c 

and  the  expressions  for  £2  and  can  be  obtained  by  cyclic 
permutation  of  the  subscripts  in  Eq.  (4.2).  We  shall  assume 
that  the  displacement  components  u^  vary  with  time  as  exp(-iio’)  , 
so  that  u ~ exp(iujt)  and  u - iouu  . Then  Eq.  (4,2)  can  be 

Of  Of  Of 

rewritten  as 

<i  ' * - 2ct(^f  + ' 

1.  2 / . / 3ul  au2  \ . ( aul  au3  \\ 

- iilDp  Ct  |^U2  ^ J + U3  \ + aXi  )j  • 


(4.3) 
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We  take  up  first  the  Poynting  vector  of  the  incident  Rayleigh 
wave , £c(i).  The  displacement  field  for  the  incident  wave  is  given  by 
Eqs.  (2.17)  and  (2.18).  We  have  taken  the  direction  of  propa- 
gation to  be  the  1-direction,  so  the  only  component  of  the 
Poynting  vector  of  interest  is  C°(i)  • The  latter  is  given 


c?d>  " |a 


2 2 
2 P0^ 


cR  1 ’U 
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C,/  e 


-2Btx3 
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Cu 


- 2 + 


+ <tt  e 
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+<tte 


4-(-f) 


•It 


-(-4)M- SI- 


MM 


(4.4) 


(4.5a) 


(4.5b) 


tt 


(“if 


((‘■4)'  *(‘4)1 


(4.5c) 


We  now  consider  the  contributions  to  the  Poynting  vector 
of  the  scattered  wave  associated  with  the  different  types  of  waves 
constituting  the  latter. 


a)  Scattering  into  longitudinal  bulk  waves 


The  scattered  displacement  field  is  given  by  Eq.  (3.39a) 


The  derivatives  with  respect  to  the  coordinate  components  are 


found  to  be 


Substituting  this  result  into  Eq.  (4.3) 


in  the  limit  x 


we  obtain 


Again  using  Eq.  (3.39a) , we  find  that  this  contribution  to  the 


Poynting  vector  can  be  written  in  the  form 


b)  Scattering  into  p-polarized  transverse  bulk  waves 


For  this  case  the  derivatives  of  the  scattered  displacement 


field  components  with  respect  to  the  coordinate  components  can 


be  written  as 
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In  the  limit  x - The  corresponding  contribution  to  the 
Poynting  vector  can  then  be  expressed  as 


rup> 


32tt2pcJx3 


(4.10) 


c)  Scattering  into  s-polarized  transverse  bulk  waves. 

In  the  limit  x -•  •,  the  expressions  for  the  derivatives 
of  the  scattered  displacement  field  components  with  respect  to 
the  coordinate  components  can  be  written  in  a form  analogous  to 
Eq.  (4.9) , namely, 


auits)(x,uj) 


dx. 


X 


U(ta> 


(x,iu) 


(4.11) 


Using  Eq.  (3.41)  we  see  that  the  Poynting  vector  from  s-polarized 
transverse  bulk  waves  takes  the  form 


where 


8tt  pc^c 


(4.12) 


u<o)(K<o)„|xo3) 

di(x0;uj) 


sincpa  cos(ktxo3cos0g) 


(4.13) 


d)  Scattering  into  other  Rayleigh  waves. 

The  scattered  displacement  field  for  this  case  is  given  by 
Eq.  (3.42).  If  we  compare  this  result  with  that  for  the  Incident 
displacement  field  given  by  Eqs.  (2.17)  and  (2.18)  and  use  the 
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Rayleigh  wave  dispersion  relation,  Eq.  (2.19),  we  see  that  the 


two  displacement  fields  differ  only  in  the  replacement 

2 .. 


A - 2 


Ct  UR 

7 7 

CR  xii 


(4.14) 


and  in  the  direction  of  propagation  parallel  to  the  surface. 

Consequently,  we  can  obtain  the  magnitude  of  the  Poynting 

vector  by  taking  that  for  the  incident  Rayleigh  wave,  Eq.  (4.4), 

and  making  the  replacement  specified  by  Eq.  (4.14). 

The  next  step  in  the  calculation  of  the  inverse  attenuation 

length  is  to  evaluate  the  energy  stored  per  unit  time  in  the 

incident  and  scattered  waves.  For  the  incident  wave,  the 

(9) 

energy  stored  per  unit  time  is  given  by'  ' 


dEo  . 

^r 


‘3  »I(1> 


(4.15) 


where  Lg  is  the  dimension  of  the  sample  parallel  to  the  surface 
and  perpendicular  to  the  direction  of  propagation.  Substituting 


Eq.  (4.4)  into  Eq.  (4.15),  we  obtain 
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(4.16) 


Turning  now  to  the  energy  stored  per  unit  time  in  the  scattered 
wave,  we  present  results  for  each  type  of  scattered  wave  separately. 
For  the  longitudinal  scattered  wave,  the  energy  stored  per  unit 
time  is  given  by 
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Using  Eq.  (4.8) , this  result  becomes 
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The  integration  over  cpg  can  be  carried  out  with  the  aid  of 
Eq.  (3.39b)  yielding 
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(4.19) 

An  expression  analogous  to  Eq.  (4.17)  can  be  applied  to  the  case 
of  p-polarized  transverse  scattered  waves.  Using  Eq.  (4.10)  for 
the  Poynting  vector  and  Eq.  (3.40b)  for  u£p^ , we  obtain  for  the 
energy  stored  per  unit  time 
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For  the  case  of  s-polarized  transverse  scattered  waves,  we 
use  Eqs.  (4.13)  and  (4.17)  and  find  for  the  energy  stored  per 
unit  time 
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Finally,  for  the  case  of  scattering  into  other  Rayleigh 
waves,  we  can  write  the  energy  stored  per  unit  time  in  the  form 
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Using  Eq.  (4.4)  and  the  replacement  specified  by  Eq.  (4.14),  we 
carry  out  the  integration  over  xg  and  obtain 
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where 


6 - 


» U + "Ut ^ itt 

RI  RTRTI  R) 


2C 


The  integration  over  cpg  can  be  accomplished  with  the  aid  of 
Eq.  (3.42b)  to  yield  the  result 


(4.24) 
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(4.25) 


We  are  now  in  a position  to  calculate  the  various  contri- 
butions to  the  Inverse  attenuation  lengths.  We  consider  separately 
the  two  cases  of  mass  defects  localized  at  the  surface  “ 0 and 
of  mass  defects  uniformly  distributed  throughout  the  crystal. 

Turning  first  to  the  case  of  defects  localized  at  the  surface, 
we  consider  a rectangular  patch  of  surface  with  dimensions  and 
Lg,  respectively,  parallel  and  perpendicular  to  the  direction  of 

propagation  of  the  incident  surface  wave  and  with  n defects  per 

s 

unit  area.  The  fraction  of  the  incident  energy  radiated  into 
longitudinal  bulk  waves  is  given  by 


dE.  / dErt 
f l " LlL2ns  dt~ / ’ “dt 


(4.26) 
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The  corresponding  contribution  to  the  inverse  attenuation  length 


is 


7TO  “ V4 


(4.27) 


For  defects  at  the  surface,  xq3  - 0,  we  find  from  Eqs.  (2.18): 
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If  we  use  these  results  and  set  xQg  equal  to  zero  in  Eqs.  (3.39c) 
and  (3.39d)  for  F1(0s,xq3)  and  (F2(9s,xo3) , we  can  simplify  Eq.  (4.19) 
for  dE^/dt  and  obtain  the  result 
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(4.30 
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In  a similar  fashion,  the  contribution  of  p-polarized 

transverse  bulk  waves  to  the  inverse  attenuation  length  can 
be  obtained  from  Eqs.  (4.20),  and  (4.26)  -(4.29)  with  the  result 

2 5 4 

j ng(Am)  uu  cR 

(/tp)  16ttP2c^€ 


The  contributions  of  s-polarized  transverse  bulk  waves  and  of 
Rayleigh  waves  to  the  inverse  attenuation  length  can  be  evaluated 
in  closed  form  using  Eqs.  (4.21),  (4.25),  and  (4.26)  -(4.29). 

The  results  are,  respectively, 
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We  now  turn  our  attention  to  the  case  where  the  mass  defects 
are  distributed  uniformly  throughout  the  crystal.  It  is  now 
necessary  to  multiply  the  energy  stored  per  unit  time  in  the 
scattered  wave  by  the  concentration  of  defects  nb  and  to  integrate 
over  xq3*  The  contribution  to  the  inverse  attenuation  length  from 
longitudinal  bulk  waves  is  then  given  by 
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dSt(xo3) 


(4.34) 


The  integral  over  x^  can  be  evaluated  analytically  if  we  ignore 

the  resonance  situation  and  set  D,  (x  :uj)  - D,(x  ;uj)  - 1.  However, 

X o o o 

the  results  are  long  and  cumbersome,  and  one  is  still  left  with 
the  integral  over  0g  to  be  done  numerically.  We  content  ourselves 
with  presenting  expressions  which  give  the  order  of  magnitude  of 
the  various  contributions: 
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V.  Discussion 

When  the  mass  defects  are  localized  at  the  surface,  our 
results  for  the  various  contributions  to  the  inverse  attenuation 
length  of  Rayleigh  waves  are  proportional  to  the  square  of  the 
change  in  mass  introduced  by  the  defect  and  to  the  fifth  power 
of  the  frequency.  In  this  respect  our  results  are  in  agreement 
with  those  of  Steg  and  Klemens^  and  of  Sakuma^5,6^  Like  Sakuma, 
we  have  also  obtained  resonance  behavior  (see  Appendix) . However, 

I 

our  resonance  behavior  is  somewhat  different  from  that  found  by 
Sakuma.  He  found  resonances  when  the  mass  defect  corresponds  to 
an  impurity  atom  lighter  than  the  atom  of  the  host  lattice  it 
replaces.  We  find  the  opposite  situation — i.e.,  resonances  when 
the  impurity  atom  is  heavier  than  the  host  lattice  atom.  We 
believe  our  result  is  more  reasonable  physically  than  that  of 
Sakuma.  A light  impurity,  under  proper  conditions,  will  lead  to 
a localized  impurity  mode  whose  frequency  lies  above  the  allowed 
band  of  frequencies  for  the  bulk  crystal.  The  impurity  mode 
frequency  therefore  lies  outside  the  range  of  Rayleigh  wave 
frequencies.  It  is  difficult  to  visualize  how  there  can  be  a 
strong  resonant  interaction  between  a Rayleigh  wave  and  an 
impurity  mode  whose  frequency  cannot  equal  that  of  the  Rayleigh 
wave.  A heavy  impurity,  on  the  other  hand,  can  lead  to  a 
resonance  mode  whose  frequency  lies  within  the  allowed  band 
for  the  bulk  crystal.  In  general,  there  will  be  no  problem 
in  finding  a Rayleigh  wave  whose  frequency  is  equal  to  that 
of  the  resonance  mode.  We  interpret  the  peaks  in  the 
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quantities  | Dx  (o ; ou)  | and  |D3(o;uu)|  to  be  at  the  resonance 

mode  frequencies  of  the  heavy  mass  defect.  The  peaks  do  not 

occur  at  the  same  frequency  because  the  impurity  mode  in  which 

the  motion  of  the  defect  is  perpendicular  to  the  surface  lies 

at  a lower  frequency  than  does  the  mode  in  which  the  motion  of  the 

defect  is  parallel  to  the  surface. 

It  is  of  interest  to  evaluate  the  order  of  magnitude  of  the 

inverse  attenuation  length.  If  we  consider  silicon  with 

p - 2.5  gm  cm  ^ , cR  “ 4.9  x 10 5 cm  sec  *,  ct  “ 5.3  x 10**  cm  sec  1, 

5 -1 

and  c^  - 9.5  x 10  cm  sec  , and  assume  that  the  defects  are 

12  -2  ■•23 

localized  at  the  surface  with  n ” 10  cm  and  Am  “10  gm, 

S 

we  obtain  for  uu  - 1010Hz  the  values  - 3.1  x 1015  cm., 

t<tp)  » 2.1  x 1015  cm,  - 2.1  x 1015  cm,  and  0.77xl015cm. 

In  these  calculations,  the  integrals  in  Eqs.  (4.30)  and  (4.31)  have 

been  taken  equal  to  unity.  We  see  that  for  the  defect  concentration 

and  frequency  considered,  the  attenuation  length  is  extremely  long. 

13 

If  we  increase  the  frequency  by  a factor  of  1000  to  10  Hz,  we 

then  find  ■ 3.1  cm,  “ 2.1  cm,  “ 2.1  cm,  and 

( R) 

l " 0.77  cm.  These  lengths  are  moderately  short,  but  the 
frequency  is  much  higher  than  those  used  in  surface  wave  devices. 

In  the  above  calculations,  the  resonances  have  been  neglected. 
Including  the  resonances  could  decrease  the  attenuation  lengths 
by  an  order  of  magnitude  or  more  at  the  resonance  frequencies. 

One  can,  of  course,  also  decrease  the  attenuation  length  by 

increasing  the  mass  change  Am  or  the  concentration  n of  defects. 

s 

In  typical  situations,  however,  it  seems  likely  that  scattering  by 

(9) 

surface  roughness'  ' will  be  more  significant  in  determining  the 
attenuation  length  than  will  scattering  by  mass  defects. 
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Appendix 

In  this  Appendix  we  evaluate  the  resonance  denominators 
Da(xo  ; u> ) («  ■ 1,2  ,3)  which  enter  the  expression  for  the  scattered 
displacement  field,  Eq.  (2.12),  and  which  are  defined  in  Eq.  (2.11). 
We  have  that 

2 


d~(x_  ; w)  - i + 


A mm 


Gcr«(xo  ’xo  JU,) 


(A . 1 ) 


We  begin  by  showing  that  Gag  (x  , x ; m)  is  diagonal  in  cr  and  P , 
and  that  G^1(x  , x ; m)  ■ G22(x  » * I m)  * G3g(x  » * I m).  For  .this 
purpose  we  first  note  that  the  Green’s  function  G^  (x  , x'  ; m)  can 

ao  .4  ^ . 

be  written  as  the  sum  of  a contribution  Gap  (x  ,x  ;m),  which  is 

the  Green's  function  for  an  infinitely  extended  medium,  and  a 

contribution  AGag(x  , x ,*<«),  which  reflects  the  semi-infinite 

nature  of  the  elastic  medium  and  ensures  that  the  surface  • 0 

is  stress-free.  The  Green's  function  G*  (x  , x'  ; uu)  diverges  as 

“P 

x -»  x;,  so  that  it  must  be  treated  separately  and  specially.  It 
is  to  this  function  that  we  turn  first. 

The  Green's  function  Ga^(x  , x'  ; m)  has  the  Fourier  expansion 

o a, r ,.)  - J , <«,.  i*3*3 ')•**.  • <*-*>  . (A.2) 


(2tt)‘ 


where  the  Fourier  coefficients  in  turn  are  expressible  in  terms 
of  simpler  functions  dag  (k||  in  | XgXg ) , which  depend  on  leu  only 
through  its  magnitude,  through  Eqs.  (2.13)  and  (2.14).  The  part 
of  dgg  (k||w  |xgXg)  which  gives  rise  to  the  Green's  function 

a ^ 4 . ' ( 9 ) 

Gag(x  , x ; uu)  has  as  its  only  nonzero  elements 


40 


« , ,,  - k||-  r 1 A3  ~ a3  1 , a3'1 

d11  (k„  im  | X3X3)  — Le  "€e  J 


2 r Ixo-xll  -ajx^-x' 


2 

2<V 


(A. 3a) 


CD  I f 

C^h  id  I ^2X3 ) * 


&22  ^11  ^ I X3X3 ) 


- ikn  , • #.  r 4 

—3  sgn  (x3  - x3)  [e 


1 -0ft ' x3  "*  x3  I 
x e 


/v  r ~“t^X3~X3l  _“t^x3"x3h 


(A. 3b) 
(A. 3c) 


2Vt 


d.3^  (k||  a)  | X3X3) 


- ikll  ___  , r -at'x3-x3'  _-“t'X3“x3' 


sgn  (x3  - X3)  J^e 


d33(k||  u)|x3x3)  - " £ _ 


(A. 3d) 


2“t* 


(A . 3e) 


The  functions  and  appearing  in  these  expressions  are  de- 


fined by 


2 J 


“t,t  ■ (k«2  - 3~  ) k»  > 
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These  definitions  can  be  combined  into  the  single  equation 


‘t,l  “ 0 


2 * 

2 (U)  + in ) ^ 

c2  ' 


(A. 5) 


where  T)  Is  a positive  infinitesimal,  and  the  branch  cut  in  the 
definition  of  the  square  root  is  along  the  negative  real  axis. 


With  the  aid  of  the  representations 


and  Eqs.  (A. 2),  (A. 3),  (2.15)  and  (2.16)  we  obtain  for  the  ele- 
Bents  of  the  Green's  function  G"fl  (x  , x ' ; uu)  : 


In  the  limit  as  x -»  x#  this  integral  diverges  at  the  upper 
limit.  In  order  to  obtain  a finite  result  in  this  limit  we  must 
iapose  a cut  off  on  the  Integral.  We  do  so  by  confining  the  inte- 
gration to  the  interior  of  a sphere  of  radius  kD,  the  Debye 
sphere.  Such  a cut  off  arises  naturally  in  a lattice  theory, 
where  kD  is  of  the  order  of  the  reciprocal  of  a lattice  parameter, 
but  it  has  to  be  Imposed  in  this  fashion  in  a continuum  theory. 

If  we  utilize  the  spherical  symmetry  of  the  problem,  we  obtain 
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where  we  have  used  the  fact  that  the  frequency  u>  has  an  infinitesi- 


mal positive  imaginary  part  (see  Eq.  (Ay 5)) 

We  now  turn  to  the  contribution  to*  Da(xQ  ; id)  from  AGag(x  , x'  ; uj) 
The  latter  is  given  by 


LGaed  a f ^ .**  • (S"  -S*)4«\lV(k,»|x3xi) 


* S a<k||>  Sve(^M)  • 
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where  Ad^v(k||iD| *3*3)  !■  that  part  of  d^^kutDlxgXg)  which  arises 
from  the  semi-infinite  nature  of  the  elastic  medium.  It  depends 
on  the  two-dimensional  wave  vector  3||  only  through  its  magnitude. 
We  now  set  x - x'  and  note  that  the  o^ly  terms  that  survive  the 
integration  in  Eq.  (A. 9)  are  those  whose  Integrands  are  even 
functions  of  k^  and  kj.  In  this  way  we  obtain  for  the  only 
nonzero  elements  of  AG^Cx  , x ; i») 

\ 

2 

AGj^ (x  , X ; id)  — j*  ® j ^Adj j (k|| u# | XgXg ) + ^ *3*3^1 

( 2tt  ) 
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( 2tt  ) 
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On  carrying  out  the  angular  integrations  in  these  expressions  we 
are  left  with  the  results  that 


AG11(x  , x ; id)  - AG22(x  , x ; id)  - 


- 4^  J dk||  k||  ^ A dtl (k||  uj  | x3x3 ) + Ad22(k||iu|x3x3)} 

kD 

if 

AG33(x  , x ; .tu)  - -r^  J dk|j  k\\  Adg3  (ky  u)  | xgx3) 


(A. 11a) 


(A. lib) 


We  have  again  introduced  a cut  off  on  the  integrals  over  ky  to 
render  them  convergent  at  the  upper  limit. 

The  expressions  for  Adaa(k|i uj  1 x3x3)  obtained  from  Appendix  A 


of  Ref.  9 are 
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In  obtaining  these  expressions  we  have  Introduced  the  notation 
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Although  the  Integrals  In  Eqs.  (A. 14)  and  (A. 15)  can  be 
evaluated  analytically  when  x3  - 0,  the  results  are  cumbersome. 

We  have  therefore  evaluated  these  integrals  numerically,  as 
functions  of  (uj/uut),  for  several  values  of  X and  k^x^.  Simpson's 
rule  was  used,  with  100  divisions  of  the  integration  interval 

_3 

(0,1).  A value  of  10  was  chosen  for  the  quantity  r\,  and  the 

integrations  were  carried  out  in  complex  arithmetic.  The 

o 1 2 2 A 

square  roots  (u  - C )2  and  (u-X  C )2  were  always  evaluated  in 

such  a way  that  their  real  parts  were  positive,  which  assured  the 

satisfaction  of  Eq.  (A. 4).  The  results  of  these  calculations 

were  combined  with  that  given  by  Eq.  (A. 8),  and  substituted  into 

Eq.  (A.  1)  to  yield  Da(xo;u>).  Results  for  | (xQ ; uu ) 1 and 

lD.,(x  :uu)l  as  functions  of  m are  presented  in  Figs.  2 and  3 
o o 

3 

for  several  different  values  of  kDx3  with  X - and  Amk^/p  - 50. 

In  each  of  these  figures,  we  see  that  there  is  a resonance  peak 
whose  frequency  at  the  maximum  increases  as  the  Impurity  is  moved 
from  the  surface  into  the  bulk.  At  the  surface,  the  peak  frequency 
is  somewhat  smaller  for  the  impurity  motion  perpendicular  to  the 
surface  then  parallel  to  the  surface.  This  indicates  that  the 


^ 1 ■">--  | — ■ 
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effective  force  constant  governing  the  impurity  motion  perpendi- 
cular to  the  surface  is  less  than  that  governing  the  motion 
parallel  to  the  surface.  A similar  effect  is  found^10^  in  the 
mean  square  displacements  of  surface  atoms. 


j 


J 


TABLE  I. 
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The  coefficient  functions  (kuj)  appearing  in 
Eq.  (3.18).  The  factor  € is  € - a o+/k2. 
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Figure  Captions 


Contour  for  the  evaluation  of  the  integral  in 
Eq.  (3.30). 

Resonance  factor  as  a function  of  frequency  for 
impurity  motion  parallel  to  the  surface. 
Resonance  factor  as  a function  of  frequency  for 
impurity  motion  perpendicular  to  the  surface. 
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